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DIGEST 

A  theory  explaining  the  break" p  of  large  liquid  drops  by  canopy 
formation  has  been  derived.  The  theory  explains  the  deformations  tiiat 
take  pJacj  during  the  disintegration  of  a  drop  in  terms  of  aerodynamic, 
hydrostatic,  and  surface-tension  pressures.  The  effect  of  liquid  viscosity 
has  been  ansumed  to  be  negligible.  Therefore,  the  theory  will  not  apply 
to  liquids  of  high  viscosity.  This  theory,  insofar  as  it  has  been  possible 
to  verify  it,  is  in  good  agreement  with  breakup  experiments  performed  with 
liquid  drops. 
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BREAKUP  OF  A  LIQUID  MASS  IN  FREE  FALL  BY  CANOPY  FORMATION: 

A  THEORETICAL  STUDY 


I.  INTRODUCTION. 


When  a  liquid  nuti  disintegrates  in  free  fall  in  r.snturbulent  air, 
three  different  modes  of  breakup  may  be  obaerved  becauae  of  varying  condi¬ 
tions  of  inatabiiity  of  the  mass.  Theae  methods  of  breakup  are  the  following: 


1.  Canopy  formation  or  "bursting- ha"  effect,"  experimentally 
described  by  Magarvey  and  Taylor*  and  Blanchard.  * 


2.  Surface  stripping,  investigated  by  Wilcox,  Piatritto,  and 

Palmer.  3 


3.  Rotation  and  oscillation,  experimentally  observed  by 
Blanchard. * 


The  first  mechanism,  canopy  formation,  refers  to  the  deformation 
of  a  liquid  mass  in  such  a  way  that  the  mass  first  becomes  flattened  on  its 
leading  side  and  subsequently  inflated  somewhat  like  a  soap  bubble  until  it 
bursts.  The  second  mechanism,  surface  stripping,  refers  to  small  amounts 
of  liquid  being  continuously  stripped  off  the  surface  of  the  man  by  aerodynamic 
forces.  The  third  mechanism,  rotation  and  oscillation,  consists  of  rotation 
of  ths  entirs  mass  until  centrifugal  force  causes  it  to  split.  Two  »maller 
manses  of  approximatsly  equal  else  usually  are  formed.  This  same  type  of 
splitting  occurs  when *Ke  maos  oscillates  into  ellipscidal  shapes  undulating 
firat  in  one  plane  and  then  in  a  plane  perpendicular  to  it. 

In  order  to  derive  a  generalised  theory  of  liquid- mass  breakup 
that  will  account  for  all  of  the  mechanisms  involved,  a  aeries  of  idealised 
situations  will  be  presented  to  make  it  possible  to  analyse  one  mechanism 
at  a  time. 

II.  ANALYSIS  OF  CANOPY^ FORMATION  l.**T-HAN!SM  IN  UQUID 
BREAKUP. 

The  most  fundamental  mechanism  of  breskup  appears  to  be  canopy 
formation.  By  fundamental  it  is  meant  that  this  mechanism  can  He  anM/sed 
independently  of  th»  others  without  imputing  many  restrictions  c.r.  the  :dealiie< 


experiment  and  without  making  many  broad  assumptions.  An  analysis  of  this 
mechanism  will  be  the  first  ctr.p  in  deriving  a  generalized  theory  of  liquic  - 
mass  breakup. 

Assume  that  a  spherical  mass  (of  radius  rQ)  of  a  nonviscous 
inelastic  liquid  of  density  p.  ,  viscosity  pj,  and  surface  tension  y,  is  at  rest 
at  time  zero  (t  =  0)  when  it  begins  to  undergo  free  fall  in  a  gravitational  field 
g  through  a  nonturbolent  still  atmosphere  of  density  pa  and  of  viscosity  of  pa- 
It  is  assumed  that  the  total  mass  of  the  liquid  remains  constant  during  free 
fall  until  the  time  when  the  canopy  bursts;  i.  e.  ,  the  amount  of  liquid  surface 
stripped^  or  evaporated  is  negligible.  In  the  case  of  a  1-gru  spherical  mass, 
the  evaporation  is  calculated  to  be  less  than  0.01%  by  Langmuir's  equation.  ^ 

As  will  be  shown  later,  this  assumption  of  constant  mass  is  only 
'  alid  for  spheres  in  a  limited  size  range.  At  time  zero  the  idealized  liquid- 
mass  profile  appears  circular  as  in  figure  1.  The  mass  is  always  described 
as  a  surface  of  revolution. 

Ill  ANALYSIS  OF  LIQUID -MASS  DEFORMATION. 

As  the  sphere  begin*  to  fall,  the  aerodynamic  pressure  on  the 
leading  hemisphere  will  begin  to  increase;  the  entire  sphere  will  begin  to  be 
flattened.  The  flattening  of  the  leading  hemisphere  will  be  more  rapid  than 
that  of  the  trailing  hemisphere.  This  can  be  visualised  by  considering  the< 
mo  ament  of  the  leading  point  as  the  aerodynamic  pressure  increases. 

A.  Development  of  Oblate  Spheroid. 


The  profile  of  the  liquid  mass,  which  is  always  described  as  a 
surface  of  revolution,  may  be  represented  on  the  original  coordinate  axes 
(figure  1),  which  are  moving  with  the  same  velocity  as  the  entire  mass 
relative  to  the  air.  As  the  asrodynamic  pressurs  increases  with  increasing 
velocity,  the  leading  point  will  move  toward  the  trailing  point  (figure  Z). 

The  distance  traveled  by  leading  point  toward  the  trailing  point  will  be 
called  "s."  As  the  leading  po.nt  move,  tcvnrd  the  trailing  point,  the  liquid 
in  the  lower  hemisphere  will  be  displaced  outward  fr-m  the  waist  of  the 
liquid  mass  along  the  path  of  leaet  resistance.  This  will  continue  to  occur 
until  the  hydrostatic  pressure  at  the  waist  of  the  mass  plus  the  surface- 
tension  pressure  caused  by  the  curvature  at  the  trailing  point  is  equal  to  the 
•  urlace  tension  at  the  waist  of  the  mass.  This  pressure  relationship  the 
iiauid  mess  may  he  expressed  by  the  following  mathematical  model: 
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ZP 

8  W 


=  Ph  + 


•t 


(i) 


where 

P  =  surface-tension  pressure  at  any  point  on  waist 

P  -  hydrostatic  pressure  at  any  point  on  waist 
n 

P  =  surface-tension  pressure  at  trailing  point 

Using  the  equations  derived  in  appendix  A,  these  parameters  may 
all  be  expressed  terms  of  s  and  rQ: 


|_L+_LI  =  v 

3(2r  -  s) 
r 

I  x*  r  I  ^ 

1  1  2 1 

8r*  +  2r  (2r  -  s)^  -  3r2  (2r  -  t*)2 

L  ©  OO  o  o 

8r4  ♦  2r  (2r  -  s)3 
o  o  o 


-  3r  {2r 


2 

•)  “  2 


3{2r  -  s) 


Ph 


f  3(2r  -  s)2  -> 

2y! — 5 — 2 - r 

I-  4r0  +  (2r0  *  •>  J 


where 

r^  =  radius  of  curvature  of  trailing  point 
r^  =  major  radius  of  curvature  at  waist 
r^  =  minor  radius  of  curvature  at  waist 
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FIGURE  2 

PROFILE  OF  LIQUID  MASS  WHEN  s  >  s  >  r 

a  o 
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The  pressure  relationship  of  the  liquid  mass,  equation  (l),m!iy 
now  be  written: 


2yIt:  +  Trj- pi  {2ro-s,g*^ 

l  2>  k 


(la) 


2Y 


3(2r  -  s)‘ 
o 


8r*  ♦  2r  (2r  -  s)3  -  3r2  (2r  -  s)2 
o  o  o  o  o 


8r*  ♦  2r  (2r  -  s)3  -  3r2  (2r  -  s) 

O  Q  O  O  O 

3(2r  -  s)2 
o 


2  -i  4 


=  Pj  (2ro 


•  )  g  +  2y 


r  3<2r«  -  »>2 

o 

4r3  ♦  (2r •  -  s)3 
^  o  o 


(lb) 


In  order  for  the  liquid  mass  to  be  stable,  an  aerodynamic  pressure 
of  T  A  must  be  equal  to  the  lesser  of  the  two  total  pressures  shown  to  be  a 
surface-tension  pressure  at  the  waist  of  P§w  and  a  hydrostatic  pressure  of 
plus  a  surface-tension  pressure  at  the  trailing  point  of  Pat.  The  aero¬ 
dynamic  pressure,  P^,  at  the  leading  point  can  be  expressed  by  the  equation 
Pa  3  l/2pav^,  where  v  is  the  velocity  of  the  air  relative  to  the  mass.  By 
solving  equation  (lb)  for  s,  the  general  appearance  of  the  liquid  mass  can  be 
determined.  The  terminal  velocity  of  the  mass  can  then  be  determined  from 
the  terminal-velocity  equation  derived  in  appendix  B.  If  the  aerodynamic 
pressure  at  terminal  velocity  is  less  than  the  hydrostatic  pressure  when 
equation  (lb)  is  satisfied,  the  mass  will  not  break  up  by  canopy  formation, 
although  it  might  disii...^rhtk  ~y  the  *">rface-btrlpping  mode. 

In  the  first  stage  of  mass  deformation,  an  oblate  spheroid  is 
formed.  When  r0,>  s  >  0,  the  radius  of  curvature,  rc,  of  the  leading  point 
will  go  to  infinity.  When  the  acceleration  at  the  leading  point,  s",  is 
expressed  in  terms  of  the  surface-tension  pressure  and  aerodynamic  pressure, 
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in  equation  is  obtained  giving  the  distance  the  leading  point  travels  ati  a  function 
of  time: 


Pi-ro  ‘ 


p  V  , 

.  +  ii . 

2  r 

c 


rl  r2 


(2) 


Rewriting  the  equation  for  v^,  derived  in  appendix  B  in  terms  of  t,  t,  and  rQ, 
the  following  equation  is  obtained: 


2  8Vo« 

'  *  3p  Cn 

a  D 


r  -  4 


.3 


dr  ♦  2r  (2r  -  >)  -  3r  (2r  -  a) 
o  o  o  o  o 


2  -i 


-1 


3<2r  -  s)‘ 
o 


tanh 


[  3pagCol 

A 

1 

8r*  +  2r  (2r  -  s)3  -  3r2  (2r  -  s)2 
o  o  o  o  o 

-  4plro3 

3(2r  -  a)2 
o 

(2a) 


The  radius  of  curvature  of  the  leading  point,  r  ,  is  expressed  in  terms  of  s 
and  rQ  in  equation  (4),  appendix  A,  as: 


rf  (2r  -  a)3  ♦  5(2r 
o  o  o 


.2  . 

,)  r«  + 
o  o 


S)3  (2r  -  S)2 
o 


3(2r  -  s)2  (r  -  s)2 
o  o 


(2b) 


Equation  (2),  which  shows  the  acceleration  on  the  leading  point, 
can  now  be  rewritten  in  terms  of  s,  x^,  and  t,  to  obtain  t  as  a  function  of 
s"  and  s. 


The  tirr»“  for  the  development  of  the  oblate  spheroid,  the  first 
steps  of  deformat:on,  can  oe  dc.ci.rur.ed  in  thie  manner.  If  equation  (2)  is 
rewritten  in  terms  of  s  and  t,  the  following  equal?  sn  is  obtained* 
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2 

\L  «Vo  J 
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P  o  o  o 

2<2'o  -  *>* 


Plt2ro  -  •) 


3(2r  -  s)2(r  -  s)2 
o  o 
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1 
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p  o 


p  p 


-2 


8r  ♦  2r  (2r  -  •)  -  3r  (2r  -  •) 

P  OP  o  o 


2^  4 
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(2c) 


By  aolving  the  above  equation  for  t  when  rQ  >  a  >  0,  the  time  to 
reach  any  atage  of  deformation  in  thia  interval  of  a  can  be  estimated. 


B.  Development  of  Disk. 


The  development  of  the  diak,  the  second  stage  of  deformation, 
occurs  when  aa  >  a  >  rQ,  where  sa  equals  the  value  <f  a  when  the  pressure 
relationship  of  the  liquid-.  * '*  equation  (lbj  is  satisiied.  In  thia  stage  of 
deformation,  the  bottom  of  the  mass  has  become  liaitoned  and  rc  *oo  .  The 
second  term  in  equation  (2)  drops  out,  and  the  following  equation  is  cbtained: 
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2p-( 2r  -  s) 

M  sj 


\t (2ro  -  *> 


(3) 


14 


Equation  (3)  can  row  be  rewritten  ir.  terms  of  s,  rQ,  and  t,  and  is 
analogous  to  the  acceleration  equation  (2:).  By  solving  these  equations  for  t, 
the  time  to  reach  any  stage  of  deformation  in  the  interval  sa  >  s  >  0  can  be 
calculated. 

C.  Development  of  Canopy. 

The  canopy  will  begin  to  form  (the  third  stage  of  deformation) 
after  «  =  sft,  as  indicated  in  figure  3.  The  minor  radius  at  the  waist,  r^, 
causes  in  effect  a  torus  of  tnickness  2  r  j.  This  will  determine  the  geometry 
of  the  canopy.  The  outer  diameter  of  the  torus  is  2  rj.  The  canopy  can  then 
be  visualised  as  an  inflating  soap  bubble  formed  from  a  soap  film  stretched 
across  a  ring.  The  ring  ir  the  case  of  the  canopy  is  the  torus.  The  body 
radius  of  the  torus,  r  j,  is  given  by  the  following  equations,  which  are  derived 
in  appendix  C: 


rT 


when  r  >  s  >  0 
o 


when  2r  >  s  >  r 
o  o 


The  inner  diameter  of  the  torus,  2rw,  is  given  by  the  following  equation: 

2rw=Mrl*rTl 


The  .nner  radius  of  the  torus,  rw,  is  constant  until  s&  «  rw  as  indicated  in 
figures  2,  3,  ard  4,  where  s&  is  the  distance  that  the  central  point  of  the 
canopy  travels  auer  s  *  s#.  When  oo  >  sa  >  rw,  then  rw  *  rc. 


t 


\ 


When  tin,  •  d»r  condition  has  been  met,  the  torus  itself  will 
expand  as  the  canopy  expands.  The  inner  xa<Uu»  of  curvature  of  the  canopy, 
rc,  can  now  b>  expn  seed  in  terms  of  rw  and  sa  u<  obtain  the  following  equation 
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2 

r 
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2s 
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I  Canopy  is  a  hemisphere 


T 


ri ?rc> 

au»avera (y  thickness  oi  canopy 
.OUKE  4 


PROFILE  OF  LIQUID  MASS  tfHEN  s  = 

a 

(Original  axes  transformed  so  that  y  = 
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The  preceding  aquation  shows  that  rc  will  be  a  minimum  when 
s  »  r.#.  Figures  2,  3,  and  4  show  the  general  appearance  of  the  liquid  mass 
when  s0  >  s  >rQ,  a  »  <>v  and  «  rw.  respectively. 

An  equation  similar  to  acceleration  equations  (2)  and  (3)  must 
now  be  found  to  det'^ibe  the  movement  of  the  central  point  of  the  canopy. 

The  acceleration  of  the  central  point  of  the  canopy  is  a  function 
of  aerodynamic  pressure  and  the  changing  radius  of  curvature.  During  the 
major  portion  of  the  interval  rw  >  »4  >  0,  the  canopy  hat  two  surfaces  with 
almost  equal  curvatures.  The  acceleration  of  the  central  point  is  now 
approximated  by  the  following  equation: 

2  2.2 

p  v  w  r  4y  *  r 

,M  s  — £ - —  .  - —  (4) 

a  2M  Mr  1  ' 

c  c  c 

where 


M  *  mass  of  tna  canopy  (cons' ant) 


a  3 

2  •'.*'1 
v  - - j - 

3p  r  C. 

a  w  D 
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'  ».«CD  V 
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'  r  (t.  4t) 
w  i 


where 


■  value  of  t  when  s  ■  s 


Substituting  the  equations  lor  Mc  and  v*  into  equation  (4),  the 
following  equation  is  obtained: 
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By  solving  equations  (2),  (3),  and  (4),  it  is  possible  to  determine 
the  time  it  takes  the  liquid  mass  to  deform  up  to  the  time  when  sa  =  rw- 


D.  Development  and  Expansion  of  Torus  and  Canopy. 

The  maximum  radius  of  the  liquid  mass,  rw,  is  no  longer  constant 
after  sa  =  r w,  but  is  equal  to  rc  ,  and  the  torus  begins  to  expand.  In  tins 
fourth  stage  of  deformation,  the  canopy  will  continue  to  inflate  until  the 
aerodynamic  fcrces  are  greater  than  the  cohesive  forces  of  the  canopy  and 
the  canopy  ruptures. 


The  volume  of  the  canopy  remains  constant  if  the  total  mass 
remains  constant.  The  canopy  will  become  increasingly  thin  as  it  undergoes 
''xpansion.  The  average  thickness  of  the  canopy,  su>  as  shown  in  figure  4, 
is  given  by  the  following  equation: 


3  volum*  oi  canopy 
u  area  of  canopy 


The  term  Ac  is  a  function  of  rc,  shown  in  appendix  D,  whero  the  following 
equations  are  derived: 
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Whan  the  conditions  of  the  following  equation,  derived  in 
appendix  £1,  are  met,  the  canopy  will  rupture  (figure  5): 


2  2 

4*r  e  C.  =  *r  p  v 
c  u  h  c  a 


(5) 


where  C,  is  the  tensile  strength  of  the  liquid. 

n 

Simplifying  equation  (5),  expressing  it  in  terms  of  sa  and  rw, 
and  substituting  Mc/p.  for  in  Vc  =  •UA<;,  the  following  equation,  derived 
in  appendix  E,  is  obtained: 
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when  .  >  s  >  0,  and 
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when  oo2  ia>  rN  and  where  rTis  taken  at  s  =  sa  and  v  is  the  velocity  of  the 
mass  just  prior  to  breaxut,. 

,  when  sa  j>  rw  £  oo ,  rw  ■  rc  and: 
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IV.  SUMMATION  OF  ANALYSIS. 


Previous  sections  of  this  report  show  that  from  the  acceleration 
equations  (2),  (3),  (4).  and  (6)  it  is  possible  to  determine  the  time  it  takes 
ior  the  mass  to  reach  any  stage  of  deformation  from  time  zero  to  the  time 
when  the  canopy  bursts. 


If  one  assumes  the  shock  of  the  liquid  disintegration  does  not 
produce  oscillations  that  would  cause  the  torus  to  break  up  (since,  at  this 
point  in  the  study,  it  appears  that  considerable  energy  is  expended  during 
breakup  of  the  canopy),  then,  after  the  canopy  has  burst,  the  torus  will 
remain  intact  until  the  aerodynamic  and  surface-tension  pressures  break  it 
up  into  smaller  more  stable  pieces.  Since  the  velocity  of  the  mass  is  known 
at  the  time  ol  ^*'ire,  it  is  possible  to  calculate  the  kinetic  energy 

of  the  torus  at  this  instant.  The  surface  energy  of  the  torus  also  may  be 
calculated.  The  torus  will  break  up  into  the  nui.iud  r*  rnece*  that  will  give 
a  minimum  of  surface  energy.  The  most  probable  number  of  pieces  wi,:,i  be 
found  when  the  torus  breaks  at  interval*  of  4.  5  diameters,  since  this  will 
give  a  minimum  of  surface  area  and,  therefore,  a  minimum  surface  energy. 
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I 

FIGURE  5 


PROFILE  OF  LIQUID  MASS  JUST  PRIOR  TO 
CANOPY  RUPTURE  WHEN  >  rw 


The  radius  of  the  spheres  into  which  the  torus  breaks  (rg)  is 


given  by  the  following  equation: 


4  3«  2  3 

3  *rs=  2  rTWrT  ;  rs=  2  rT 


where 


2  3  volume  of  torus 


iir  r_  -  r_  - 
c  T  T 


=  0  and  r^,  is  taken  just 


before  disintegration. 


If  r8  is  less  than  the  maximum  stable  drop  size  determined  by  the 
acceleration  equations  (2),  (3),  (4),  and  (6),  the  system  has  become  stable. 

If  rg  is  more  than  the  maximum  stable  size,  the  whole  process  outlined  in 
this  paper  must  be  repeated  for  each  drop,  taking  into  consideration  that  the 
initial  velocity  of  the  air  is  no  longer  zero. 

V.  DISCUSSION. 


When  deriving  this  theoretical  treatment  of  canopy -formation 
b~eakup,  two  assumptions  were  made.  First,  it  was  assumed  that  the  mass 
of  the  liquid  remained  constant  throughout  the  deformation  process.  This 
assumption  is  valid  only  for  relatively  small  liquid  spheres.  It  follows  from 
the  definition  of  surface-tension  pressure  that,  as  a  liquid  mass  becomes 
distorted  into  a  larger  shape,  the  surface-tension  forces  holding  the  outer 
layers  of  liquid  will  decrease  with  increasing  surface  area.  The  time 
required  to  reach  any  given  free-fall  velocity  and  the  related  aerodynamic 
force  will  also  decrease  as  the  sphere  be  omes  larger.  Consequently,  at  • 
any  given  instant  in  the  dsformation. process,  there  will  be  a  greater  prob¬ 
ability  of  stripping  off  of  surface  layers  of  liquid  from  the  mass.  For  - 
spherical  liquid  masses  in  the  diameter 
a  negligible  amount  of  surface  stripping, 
surface  stripping  longer  be  neglected,  3  and  the  theory  derived  in  this 

paper  must  be  combined  with  a  tneory  of  liquid-mass  breakup  caused  by 
surface  stripping.  A  future  study  is  contemplai.d  to  «how  tha*  very  large 
liquid  masses  break  up  almost  entirely  by  surface- stripping  effects. 

The  second  assumption  was  that  the  effect  of  the  viscosity  of  the 
liquid  was  negligible.  The  viscosity  may  be  neglected  in  the  cauc  of  so*»ll 


range  of  1. 2  cm  to  4.  0  cm,  there  is 
*  *  For  diameters  of  over  4.  0  cm, 
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masses  of  relatively  nonviscous  liquids  of  high  surface  tension  such  as  water, 
but  cannot  be  neglected  in  general.  The  effect  of  increased  viscosity  will  be 
to  slow  down  the  rate  of  deiormation  by  decreasing  the  acceleration  of  the 
leading  point  and,  possibly, to  slow  down  the  acceleration  to  the  point  where 
the  mass  will  be  stable;  i.  e.  ,  the  mass  will  have  reached  terminal  velocity 
before  the  aerodynamic  force  can  overcome  the  viscosity  and  surface  tension 
of  the  liquid.  Consequently,  the  rate  at  which  a  liquid  mass  will  break  up  can 
be  decreased  by  increasing  the  viscosity,  or  it  can  be  increased  by  decreasing 
the  surface  tension.  As  a  liquid  mass  becomes  sufficiently  large  so  that 
viscosity  must  be  considered,  the  main  mechanism  of  breakup  is  no  longer 
canopy  formation  but  surface  stripping. 

This  theory  was  checked  and  found  to  be  in  agreement  with  the 
studies  reported  by  Magarvey  and  Taylor*  and  Blanchard.  2  These  agreements 
are  shown  in  appendixes  F  and  G.  The  difficulty  of  solving  acceleration 
equations  (2),  (3),  (4),  and  (6)  analytically  has  made  it  impractical  to  check 
the  times  of  deformation.  The  predicted  general  geometry  assumed  by  the 
liquid  mass  is,  however,  in  very  good  agreement  wit!  experimental  results. 

The  coefficient  of  cohesion  between  the  liquid  surfaces  is  uncertain  since  it 
depends  on  the  impurities  in  the  liquid.  The  coefficient  of  cohesion  can  be 
calculated  by  combining  this  theory  with  data  obtained  in  canopy-breakup 
experiments.  This  cohesion  coefficient  is  shown  in  appendix  G  to  be  consid¬ 
erably  less  than  1  atmosphere,  which  is  much  lower  than  would  be  predicted 
from  molecular  considerations.  The  effect  of  impurities  on  the  tensile  strength 
of  a  '  iquid  is  still  not  known  in  a  quantitative  sense.  Before  this  theory  can  be 
tested  fully  or  put  to  some  practical  use,  more  work  on  the  tensile  strengh  of 
liquids  is  necessary.  A  satisfactory  method  for  solving  nonlinear  differential 
equations  of  the  second  order,  such  as  those  describing  the  acceleration  of  the 
trailing  point  of  the  liquid  mass,  must  also  be  found.  For  the  time  being,  a 
numerical  approximation  must  suffice. 

In  conclusion,  it  should  be  mentioned  that  the  mathematical  model 
postulated  for  the  deformation  process  is  certainly  not  an  exact  one.  It 
appears  to  be  a  very  <-h.«e  approximation  as  shown  by  figures  and  examples  in 
appendixes  C  and  F  V  ellipsoids  oi  revolution  had  been  postulated  as 
mathematical  models  instead  of  sections  of  spheres,  theoiy  would  not 
have  been  in  as  good  agreement  with  experimental  results  as  it  is. 
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GLOSSARY 


A  -  maximum  cross-sectional  area  of  liquid  maaa 

-  area  of  in-  *r  aurface  of  canopy 

Cp  -  coefficient  of  drag 

r  -  maximum  aerodynamic  force  acting  on  liquii  maaa 

max 

g  -  acceleration  caused  by  gravity 

M  -  total  maaa  'f  liquid 

-  maaa  of  canopy 

-  aerodynamic  preaaure  on  leading  point  of  liquid  maaa 

-  hydroatatic  preaaure  at  any  point  or  at  wait l  of  liquid  maaa 

P>w  -  aurface -tenaion  preaaure  at  any  point  on  waiat  of  liquid  maaa 
Pat  -  surface -tenaion  preaaure  at  trailing  point 

Pq.  -  preaaure  of  atability 

rQ  -  original  radiua  of  spherical  liquid  mass  (figure  1) 

r  j  -  major  radius  of  curvature  at  waist  of  maaa  (figure  2) 

r £  -  minor  radius  of  curvature  at  waist  of  maaa 

r  -  radiua  of  curvature  c'  ceni.i.'  leading  point  of  liquid  mass 

c 

-  radiua  of  curvature  of  central  trailing  point  of  liquid  maaa 

r^.  -  one-half  thickness  of  torus;  i.  e. ,  radius  of  body  of  torus 

r  -  inner  radius  of  torus 

* 
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r#  -  radius  of  spheres  into  which  torus  breaks 
R^  -  Reynolds  number 

s  -  distance  leading  point  travels  toward  trailing  point  in  moving 

coordinate  system,  where  trailing  point  moves  with  same  velocity 
an  coordinate  axes  (figures  2  and  3) 

-  value  of  s  when  equation  (lb)  is  satisfied 

s  -  distance  central  point  on  leading  surface  of  canopy  travels  after 

equation  (lb)  is  satisfied 

s^  -  average  thickness  ot  canopy  (figure  3) 

t  -  time  elapsed  between  any  given  intervals  of  s  or  e^ 

t,  -  value  of  t  when  s  =  s 

i  a 

t..  -  total  time  elapsed  from  s  >  0  to  s  *  r 

u  aw 

v  -  velocity  of  atmosphere  relative  to  liquid  mass 

v  -  terminal  velocity  of  mass 

Vj  -  volume  of  liquid  above  waist  of  mass 

V  -  volume  of  liquid  below  waist  of  mass 

£ 

Vc  -  volume  of  canopy 

p^  -  density  of  air 

pj  -  density  of  liquid 

-  viscosity  of  air 

Pj  •  viscosity  of  liquid 

% 

y  -  surface  tension  of  liquid 
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APPENDIX  A 


DERIVATION  OF  RADI!  OF  CURVATURE  AT  WAIST  AND 
AT  LEADING  AND  TRAILING  CENTRAL  POINTS 


In  the  rungs  rQ  >  sj>  0,  where  rQ  ia  the  original  radius  of  the 
spherical  mass  and  s  is  the  distance  that  the  central  leading  point  travels 
toward  the  central  trailing  point  ip  a  moving  coordinate  rystem  where  the 
central  trailing  point  remains  stationary  relative  to  the  coordinate  axes, 
there  exists  the  following  relationship: 


4 

3 


it  r 


V  4  V 
12 


where 

Vj  -  volume  of  liquid  above  waist  of  mass 

V  =  volume  of  liauid  below  waist  of  mass 

Vj  and  V^,  described  sections  of  spheres,  are  volumes  of 
revolution  and  may  be  visualised  in  figure  1 . 


The  volumes  V^  and  V^  then  are  calculated  readily  in  terms  of 
the  radius  of  curvature,  r^>  of  the  central  trailing  point  and  the  radius  of 
curvature,  rc,  of  the  central  leading  point. 


In  general: 


When  r  >  s  >  0: 
o 


V-**'Vo 


(1) 
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In  general: 
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3. 


Since  the  mu»  of  the  liquid  ia  aaaumed  conatant  and 


—  ir^  =  Vj  4  V2>  the  following  relationahip  exiata: 
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M  •4,0*  •  <2ro  -  »>3  ro  *  5,o  <2r„  -  ‘I2  ♦  <»•„  -  «>3  1»«„  -  ‘i2  4 

-  ->2 1'0  -  -I2  Hl 


r,  *  maximum  radius  of  the  mass  (figure  1) 


By  expressing  in  terms  of  its  position  in  the  moving  coordinate  system, 
the  following  relationships  are  obtained: 


2  2  2 
x  +  y  a  r. 


2  2 

*  .rj 


y  a  r.  -  r 
k  o 


.22.  ,2  ,  2 

•  •rl*rk-{rk-ro)  =2Vo*ro 


r .  *  2rr  -  r 
1  k  o  o 


for  r}  is  derived: 


After  substituting  equation  (3)  for  r  ,  the  following  expression 

_ _  i i.  * 


_  j  •  ‘>3  -  ^  -  -I2 12 

L  512'.  -  'I2 


r£  »  minor  radius  of  curvature  at  the  waist  of  the  mass 

When  r.  ■  r  ,  r,  »  r  and  as  r,  — >  ,  r ,  — »  0. 

t  o  t  o  1  2 

Therefore,  by  symmetry,  the  following  expression  for  r_  is  derived  (figure  2): 
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APPENDIX  B 


DERIVATION  OF  VELOCITY  OF  AIK  RELATIVE  TO  FREE-FALLDIG  MASS 


7‘  e  total  acceleration  acting  on  a  mass  may  be  expressed  by  the 
following  equation; 


dv 

dt 


=  g 


p  v  A  C_ 

a _ D_ 

2  M 


where 


p^  =  density  of  air 
g  =  acceleration  due  to  gravity 

v  =  velocity  of  air  relative  to  mass 

„  2  . 

A  sir  «  maximum  cross-sectional  area 

t  *  time  elapsed  after  beginning  of  free  fall 

C0  =  coefficient  of  drag  of  mass  (assume  constant) 

4  3 

M  =  Pj  y  w  rQ  3  total  mass  of  sphere  (assume  constant) 

By  solving  equation  (1)  for  v,  the  following  expression  is  obtained: 
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By  rewriting  and  squaring  equation  (2),  the  following  expression  is  obtained: 


v2  =  a2  tanh2(at/k) 


(3) 


By  substituting  the  original  values  for  a  and  k  in  equation  (2), 
expression  for  the  .  quare  of  the  velocity  is  obtained: 
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8p 


3 
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tanh 
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the  following 


(4) 


where 

Sr*  ♦  2r  (2r  -  a)3  -  3r2J2r  -  s)2 

O  O  O  O  O 

3(2r  -  a)2 
o 

r^  =  original  radius  of  spherical  mass 

s  *  distance  traveled  by  leading  point  toward  trailing  point 
Pj  =  density  of  liquid 


The  coefficient  of  drag,  C^,  is  not  constant  for  a  sphere,  but 
dec  eases  with  increasing  Reynolds  number  at  low  Reynolds  numbers.  The 
mass,  however,  does  not  remain  spherical.  It  assumes  the  shape  of  a  section 
of  a  sphere  and  then  becomes  parachute -shaped.  It  is,  therefore,  convenient 
to  assume  a  constant  coefficient  of  drag  of  about  0.  4,  which  should  give  a 
fairly  good  approximation  of  the  terminal  velocity  of  masses  who*e  Reynolds 
number  is  about  10,000  at  terminal  velocity.  The  assumption  of  a  drag 
coefficient  of  about  0.  4  is  justified  and  can  be  seen  by  checking  a  graph  of 
drag  coefficients  versus  the  Reynolds  numbers  of  spheres,  cylinders, 
hemispheres,  and  disks:  shapes  that  approximate  the  geometry  of  the  deform¬ 
ing  liquid  mass.  The  g.ua.cst  '.rag  .  ffirient  error  will  be  encountered 
when  the  Reynolds  number  is  less  than  100.  (This  er~nr  is  very  small  at 
Reynolds  number  greater  than  100.)  Consequently,  the  equation  derived 
above  should  yield  reasonable  values  for  masses  that  have  a  Reynolds  number 
of  less  than  100  for  a  relatively  short  period  of  time  (less  than  10%)  during 
free  tall. 
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APPENDIX  C 


DERIVATION  OF  BODY  RADIUS  OF  TORUS 


’  figure  1,  r-p  ia  the  radiua  of  the  body  of  the  torua  and  r 2  the 
effective  radiua  of  curvature  at  the  waist  of  the  mass. 

By  trigonometric  analogy,  when  r  AaiO 


in  tigure  2. 

in  figure  3. 


-IH 

\  1 1 


When  a  >  rQ,  the  liquid-maaa  profile  haa  the  appearance  shown 


Figure  2  may  be  approximated  by  the  mathematical  model  shown 


By  trigonometric  analogy  r  j  may  be  expressed  as: 


r  _  * 

T 


2r  -  a 
o 


when  2r  >  a  >  r 
o  o 


FIGURE  2 


ACTUAL  APPEARANCE  OF  LIQUID  MASS  WHEN  2r  >  s  >  r 

o  o 
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V  / 

- 

H — 2r2 


FIGURE  3 


APPROXIMATE  APPEARANCE  OF  FIGURE  2 
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APPENDIX  D 


DETERMINATION  OF  CANOPY  AREA 


To  determine  the  area  of  the  canopy,  the  profile  of  the  canopy 
ia  given  in  figure  1  aa  a  section  of  a  circle  of  radius  r  . 

2  2 
r  ♦  s 

W  a 

r  -  — r  1  ■  ,  where  r  *  r,  -  r_  when  s  20  (i  is  the  distanci 

c  2s  w  1  T  a  a 

a 

traveled  by  the  central  point  of  the  canopy  in  the  moving  coordinate  system; 

s  =  y  after  s  =  s  ). 
a  a 

The  area  of  the  canopy,  A  ,  is  obtained  by  rotating  a  section  of 
the  circle  around  the  y-axis. 

dA^  =  2  ir  x  ds,  da  =v4x^  -  dy^ 


2  2  2 
x  +(y-s  ♦  r  )  -  r 
a  c  c 


(equation  of  a  circlt 


2  2  2 

..  x  ♦  y  *  2s  r  -  s  -  2y(r  -  s  ) 
a  c  a  c  a 

,  ,  2.2 

,2.2  ,  r  dx 

,  ,  2  x  dx  ....  2  c 

and  -dy  *  -  2  •  *°  that  “*  =  — j.'  2 

r  -  x  r  -  x 

c  c 

.  .  When  r  >  s  >  0 

w  a 


a  ■> .  ;  w  2  2.  2  , 

A  »  2»r  J  >«•-  rr  )  dx 

c  C  0  c 


1  lrw 


2  2  2 
A  *  2vr  -  (r  -  x  )  ‘ 
c  c  c 


When  s  >  r 
a  w 


• 

2  r'c  ,  2  2  ~1  , 

A  *  2«r  ♦  2»r  x(r  -  x  )  dx 

c  c  c  J  c 

r 

w  i 

2  2  2  2 

A  =  2irr  +  2wr  (r  -  r  ) 
c  c  c  c  w 

Ac  =  2"c  [rc  *  (rc  •  rt>  2  ] 


1 


A  ie  the  ere*  of  the  inner  surface  of  the  canopy.  The  area  of 
both  surfacesCis  approximately  2  A^.  As  >  this  approximation 

becomes  exact. 
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APPENDIX  E 


DETERMINATION  OF  DISINTEGRATION  POINT 


The  canopy  will  rupture  when  the  aerodynamic  force  becomes 
greater  than  the  tenaile-atrength  force  holding  the  canopy  intact.  The  point 

given  by  the  following  expression: 

when  s  £  r  (1) 

aw  ' 

C,  *  tensile  strength  of  liquid 
h 

=  density  of  atmosphere 
v  a  velocity  of  atmosphere  relative  to  maes 

Simplifying  and  solving  equation  (1)  for  r  ,  the  following  equation 

is  obtained: 


of  equilibrium  between  the  two  forces  is 

2 

,  ,  2  V 

2*r  s  C.  *  »r  *  . . 

c  u  h  c  2 


where 


r^  *  radius  of  curvature  of  canopy 


s  a  average  thickness  of  canopy 


r  » 
c 


4s  ck 

u  h 


V 


2,  2 

s  ♦  r 

a  w 

Since  r  »  —  — v  ■  -  ■ .  it  is  possible  to  find  the  value  of  s  at 
c  * 

a 

which  the  canopy  will  rupture  if  s^  is  also  expressed  as  a  function  of  s^. 


Since  a  s  — —  a  .volume  of  canopy 
u  Ac  area  of  canopy 


4V  ch  4V  ch 

At-  •  -c  ■■ —  ,  or  A  r  =  L,  letting  — — — ; —  =  L. 
c  c  2  c  c  £■ 

p  v  p  v 

a 


From  append**  D,  when  ■  >  r  : 
rr  aw 


r  A  *  2ir 
c  c 


r  2  2  - 

r  4  s 
w  a 

<• 

r  2  A  2  *,224 

w  a  .  /  w  w  a  a  2  ] 

2  i 

2« 

2s  +  .2  *  rw 

a 

a  '  4«  / 

L  a 

.  \  r  A  =  f(s  ) 
c  c  a 


obtained: 


By  solving  equation  (1)  in  terms  of  r^,  the  following  equation  is 


4  2Lrc  *  L2  „ 

rc - 2—  ♦  -j-  =  0 


(2) 


If  equation  (2)  is  solved  for  r  and  then  substituted  in  the  equation  for  sa, 
the  point  of  canopy  rupture  isC obtained.  The  velocity,  v,  can  be  assumed 
constant  at  this  point  and  is  approximated  by  the  following  equation: 


v  * 


I 

'  3p  f2C 

a  w 


1 

2 
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APPENDIX  F 


SAMPLE  CALCULATIONS  OF  CRITERION  FOR  STABILITY  OF  DROPS 


ConsiCwr  a /spherical  drop  of  water  for  which  rQ  =  0.  45  cm, 

.y  *  surface  tension  =  72  dynes/cm,  and  p|  =  1  gm/cc.  Solving  the  acceler¬ 
ation  equation  (lb)  of  the  t4xt  for  sa  and  then  substituting  sft  into  the  appropriate 
equations  derived  in  appendixes  A  and  C,  the  following  values  are  obtained. 


s  *  0.  435  cm 
a 

r  =  0.  656  cm 


(2r  -  s  )  =  0.  465  cm 

O  ft 

r_  =  0.  216  cm 
T 


r^  *  0.  719  cm 

r  *  0.  440  cm 
w 


2\ 


-U-iJ.p, 

rl  r2J  1 


(2r  -  s  )g 

O  ft 


♦  ix 

rk 


(1) 


or  144 


1  .  0.656  1 


All  the  spherical  masses  being  considered  have  Reynolds  numbers 
in  the  neighborhood  of  10.000.  Therefore,  insofar  as  a  spherical  shape  is 
cc'  earned,  Cq  •  0. 4  is  a  good  general  value  for  the  coefficient  of  drag.  By 
using  the  equation  derived  in  appendix  B,  the  terminal  velocity  for  a  spherical 
mass  may  be  determined.  This  velocity  is  not  the  terminal  velocity  that  the 
mass  actually  reaches,  but  is  the  velocity  that  the  mass  would  reach  if  it 
retained  its  original  spherical  shape.  The  criterion  for  stability  is  that  the 
hydrostatic  pressure,  P^,  at  the  leading  point  of  the  spherical  mass  be 
greater  than  the  aerodynamic  pressure,  Pa>  of  the  spherical  mass  at  terminal 
velocity.  This  relationship  may  be  expressed  in  the  following  manner: 


PJ(2r{>)  g  > 


(««**b»Uty  criterion) 
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where 

r  =  radius  of  sphere 
o 

Pj  -  density  of  liquid 

Cq  =  0.4  drag  coefficient  (whenB^2  10,000) 

The  stability  criterion  has  the  following  appearance  when  r  =0.4! 
Pj  =  1:  p^(2rQ)g  =  883  dynes/cm^  - 

4ro2pig  ,  2  _ 

- — -  =  756  dynes/cm  =  P 

A  9- mm-diameter  mass  of  water  is  then  seen  to  be  stable  and  will 
not  disintegrate  by  canopy  formation  in  agreement  with  the  finding  of 
Blanchard^ and  Magarvey  and  Taylor.  *  This  does  not,  however,  preclude 
disintegration  by  rotation  or  oscillation.  g 

Magrrvey  and  Taylor  experimentally  determined  that  the  minimun 
waterdrop  size  that  would  disintegrate  by  canopy  formation  was  about  1. 2  cm 
diameter.  Applying  the  stability  criterion  to  a  1.  2-cm  drop,  the  following 
relationships  are  obtained: 

2 

Pj2rQg  *  1*175  dynes/cm  = 


4Vig  ,  .  2  _ 

— —  *  1, 175  dyne*/ cm  *P 

3Cd  a 


This  :s  in  excellent  agreement  with  Magarvey  and  Taylor. 
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APPENDIX  G 


SAMPLE  CALCULATION  OF  SIGNIFICANT 
PARAMETERS  FOR  WATERDROPS 


Using  the  formulas  derived  in  the  previous  appendixes  and  in  the 
body  of  the  report  and  applying  them  to  a  waterdrop  1 5  mm  in  diameter  for 
which  rQ  =  0.  75  cm,  the  following  values  for  the  significant  parametero  are 
obtained: 


2r  -  s  *  0. 490  cm 

o  a 

r,  =2.5  cm 

k 

rj  *  1 .  79  cm 

r^  -  0.  314  cm 

r *  0.  086  cm 
T 

r  »  i.  -  r_  »  1. 704  cm 
w  1  T 


equation  (lb) 
appendix  A 
appendix  A 
appendix  A 
appendix  C 


The  torus  will  begin  to  expand  when  704  cm. 


When  s  »  6  cm,  r 

a  c 


2X  2 
s  ♦  r 

a  w 


«  3.  242  cm 


When  r  «  3.  242  cm,  v 

c  max 


a  3  * 

i  8r  gp.t 

»  I - -  "’V  ]  *  298  C 

*  3p  r  / 

r»  n 


m/ sac 


Maximum  aerodynamic  force  (F  )  is  given  by  tb*  following  expression: 

max 

F  «  4  p  V*  *  r2  a  1. 8  X  to3  dynes 
a  2  a  max  c  ' 

m  A  v 


T,iis  i-  •''<*  tote1,  force  available  for  disrupting  the  canopy. 


*  *  4 


From  appendix  F,  F^  =  ^,rrc8u^i 


max 


C,  = 


h  2irr  » 
c  a 


max  ,  ,  2 

dyne./cm 


volume  of  canopy  _  c 

»  =  - m -  -  x 

c 


area  of  canopy 


4  3.2  2 

V  =  —  irr  =  aw  r  r 
c  3  o  w  T 


V  =  1.  77  cc  -  0.  242  cc  =  1.  528  cc 
c 

1 

„  2  2 .  2  1 
A  =  2ffr  r  +  (r  -  r  ) 

c  c  c  c  w  | 


A  *  122  aq  cm 
c 
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1.  528 


1. 22  X  10 


*  1. 25  X 10  ^  cm 


c  =  - _ -  *  7.  06  XI 03  dyne./. q  cm 

[6.28  (3.242)]  (0.0125)  .q  cm 


Thi.  value  of  Ch  1.  .ten  to  be  much  lower  than  would  commonly 
be  oredicted  by  quantum  mechanic,  for  water.  It  might  be  explained  that  a 
ti  i  in  the  canopy  begin,  around  a  nucleu.  that  could  he  an  impurity  in  the 
liquid. 

The  point  of  bi  eakup  wa.  cho.en  from  the  photograph,  in  Magarve 
and  Taylor. 1  Th"  1  .  *kal»*»d  parameter,  are  in  very  good  agreement 
with  the  photograph*  reproduced  in  the  ..  ove-mentioned  article. 
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